Naval  Research  Laboratory 

Washington,  DC  20375-5320 


NRL/MR/5341--98-8198 


Detection  Performance  of  Signals  in 
Dependent  Noise  From  a  Gaussian 
Mixture  Uncertainty  Class 


K.R.  Gerlach 

Target  Characteristics  Branch 
Radar  Division 


August  31,  1998 


DTIC  QUALITY  HCEPECT21D  I 


Approved  for  public  release;  distribution  is  unlimited. 


19980904  060 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMBNo.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA 22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  {Leave  BlanK) 

2.  REPORT  DATE 

3.  REPORT  TYPE  AND  DATES  COVERED 

August  31,  1998 

Interim 

4.  TITLE  AND  SUBTITLE 

5.  FUNDING  NUMBERS 

Detection  Performance  of  Signals  in  Dependent  Noise  From  a  Gaussian  Mixture  Uncertainty  Class 

PE-61 153N 

PR-LR021-05-RH 

6.  AUTHOR(S) 

K.R.  Gerlach 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

Naval  Research  Laboratory 
Washington,  DC  20375-5320 

NRL/MR/5341— 98-8198 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

800  North  Quincy  Street 

Arlington,  VA  22217-5660 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

1 2a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 

1 3.  ABSTRACT  ( Maximum  200  words) 

Detection  performance  in  terms  of  probabilities  of  detection  and  false  alarm  ( PD  and  PF)  is  derived  for  a  desired  signal  embedded 
in  dependent  (correlated)  multivariate  noise  from  a  Gaussian  mixture  uncertainty  class.  This  uncertainty  class  is  defined  using  upper 
and  lower  bounding  functions  on  the  univariate  Gaussian  mixing  distribution  function.  For  a  given  detector  operating  point,  expres¬ 
sions  for  PD  and  PF  are  derived.  In  addition,  nontrivial  lower  and  upper  bounds  of  performance  are  derived  which  are  functions  of  the 
bounding  functions  of  the  uncertainty  class,  the  detection  threshold,  and  the  desired  signal  strength  (only  for  PD). 

14.  SUBJECT  TERMS 

15.  NUMBER  OF  PAGES 

Detection 

Dependent  Noise 

19 

Signal  Processing 

Clutter 

Non-Gaussian  Noise 

16  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

UL 

NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std  239-18 
298-102 


CONTENTS 


I.  INTRODUCTION  .  1 

II.  PRELIMINARIES  .  2 

A.  Gaussian  Mixture  Detectors .  2 

B.  Problem  Statement .  4 

III.  PD  AND  PF  DERIVATION  .  5 

A.  Equivalent  GLRT  .  5 

B.  Derivation .  6 

IV.  PERFORMANCE  BOUNDS .  8 

A.  Useful  Lemmas .  8 

B.  Derivation .  10 

C.  Performance  Spread .  13 

V.  SUMMARY .  14 

VI.  REFERENCES  .  15 


iii 


DETECTION  PERFORMANCE  OF  SIGNALS  IN  DEPENDENT  NOISE 
FROM  A  GAUSSIAN  MIXTURE  UNCERTAINTY  CLASS 


I.  INTRODUCTION 

The  discrete-time  detection  of  signals  in  non-Gaussian,  dependent  (correlated)  noise  is  a  difficult 
problem  of  fundamental  interest.  Aside  from  the  non-Gaussian  aspect,  the  problem  is  complicated  by  the 
correlation  in  the  noise  data.  To  attack  the  detection  problem  requires  specification  of  a  multivariate  non- 
Gaussian  distribution  that  accounts  for  both  the  correlation  function  and  the  non-Gaussian  marginal 
distributions.  It  is  desirable  to  find  a  multivariate  distribution  of  sufficient  generality  to  model  the  desired 
noise  characteristics  (i.e.  correlation  and  non-Gaussian  marginals),  yet  of  sufficient  tractability  to  yield 
a  solution  to  the  detection  problem.  A  general  class  of  distributions  that  allows  one  to  accomplish  this 
is  the  Gaussian  mixture  distributions  (GMD). 

The  GMD  has  found  utility  in  modeling  radar  clutter  [1-6,18],  A  Rayleigh  mixture  distribution 
(RMD),  which  is  closely  associated  with  the  GMD,  is  often  used  to  model  the  amplitude  characteristics 
of  radar  clutter  [1-6].  Such  common  amplitude  clutter  models  as  Weibull,  K,  and  contaminated  Rayleigh 
can  be  represented  as  RMD’s  [1,6].  These  models  were  formulated  to  account  for  the  target-like  spikes 
that  were  observed  in  experimental  clutter  data  [2,3,7-12].  For  multivariate  clutter  modeling,  the  RMD 
has  a  natural  extension  to  modeling  clutter  as  a  multivariate  GMD  where  often  the  only  knowledge  one 
has  of  the  inputs  is  the  univariate  amplitude  distribution  and  the  correlation  matrix  [1,6,18], 

Hence  much  of  the  information  necessary  to  define  a  detection  scheme  for  a  desired  signal 
embedded  in  correlated  clutter  modeled  as  a  GMD  is  known.  For  our  purposes,  this  detection  structure 
takes  the  form  of  the  Generalized  Likelihood  Ratio  Test  (GLRT)  [13]  whereby  the  unknown  complex 
amplitude  of  the  desired  signal  is  estimated. 

Quite  often  the  univariate  mixing  distribution  of  the  RMD  (or  GMD)  is  unknown  but  can  be 
lower  and  upper  bounded  (e.g.,  using  histograms  constructed  from  measured  input  data).  Hence,  if  one 
designs  a  detector  for  a  given  GMD  which  is  within  this  class  of  bounded  distributions,  the  detection 
performance  will  vary  from  the  expected  nominal  performance  when  the  inputs  are  not  this  same  GMD 
but  still  a  member  of  this  class. 

A  minimax  detector  approach  (see.  e.g.,  [14,15])  for  finding  a  good  detector  for  a  given  input 
uncertainty  class  is  as  follows.  Find  a  distribution,  Pop ,  that  is  a  member  of  the  given  input  uncertainty 
class  on  which  we  base  our  detection  scheme  (this  is  sometimes  called  the  detection  operating  point) 

which  has  the  following  saddle  point  property;  R(Pgp,Pop)  <  R(Pop,Pop)  ^  RiP^P^).  For  these 
inequalities  /?(•,•)  is  some  measure  of  the  detector’s  performance  where  the  first  argument  indicates  the 
input  distribution  and  the  second  argument  indicates  the  detector  operating  point  distribution.  The  larger 
R  is  the  better  we  judge  the  detector’s  performance.  Pop  and  P,  are  any  distributions  from  the  input 
uncertainty  class. 

We  do  not  use  the  minimax  detector  approach  in  that  we  choose  any  distribution  from  within  the 
uncertainty  class  as  a  detector  operating  point.  Thereafter  we  quantify  the  performance  of  the  detector 
when  the  input  distribution  is  any  from  the  uncertainty  class.  The  performance  will  degrade  when  the 
input  distribution  is  not  matched  with  the  detector’s  operating  point  distribution.  Specifically,  in  this 
paper,  we  derive  the  GLRT  and  performance  measures  (probabilities  of  detection  and  false  alarm) 
associated  with  a  given  GMD  (the  designed  operating  point)  for  inputs  from  a  given  GMD  uncertainty 
class  of  bounded  distributions.  In  addition,  nontrivial  upper  and  lower  bounds  of  performance  are  derived 
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which  are  functions  of  the  bounding  functions  of  the  uncertainty  class,  the  detection  threshold,  and  for 
the  detection  probability,  the  desired  signal  strength. 

II.  PRELIMINARIES 

A.  Gaussian  Mixture  Detectors 

Consider  the  following  binary  hypotheses 


where  z,  x,  and  s0  are  AMength  vectors  belonging  to  the  A'-dimensional  complex  space,  C*,  and 
z  =  {zn;  n  =  1,  2,  ...,N}  is  a  sequence  of  observations  or  samples.  The  H0  hypothesis  corresponds  to 
the  case  when  no  desired  signal  is  present.  The  //,  hypothesis  corresponds  to  the  case  where  the  desired 
signal  vector,  as0,  is  corrupted  by  an  additive  zero-mean  noise  vector,  x.  Here  a  is  the  unknown  complex 
amplitude  of  a  known  signal  vector,  s0,  and  Sq  is  normalized  so  that  s0H  s0  =  1  where  H  denotes  the 
conjugate  transpose  operation. 

If  the  noise  is  modelled  as  a  GMD  then  z  has  the  following  TV-variate  probability  density  function 
(pdf)  under  //„  (i  =  0,1): 


f,  (z) 


1 


irN  dec 


1, 7 e,ip 


T 


dP 


(2) 


where 


<7,  = 


zH  <f>''  z  ,  i  =  0, 

(z  -  oSq)"  <t>'  (z  -  asj  ,  /  =  1, 


(3) 


$  =  normalized  covariance  matrix  of  x  (assumed  nonsingular), 

r  =  variance  of  any  element  of  x, 

P T  =  Gaussian  mixing  distribution, 
fir  =  support  of  r,  Cr  £  [0,  <»), 

and  dec  denotes  the  matrix  determinant.  The  normalized  covariance  matrix  is  defined  as  follows.  Set 
x  =  (x,,  jc2,  .  .  .  ,  xN)J  where  T  denote  the  transpose  operation.  Then 


_  £{xx»} 

E{t) 

where  E  denotes  the  expected  value.  For  the  case  where  the  N  inputs  are  equi-spaced  time  samples  from 
a  wide-sense  stationary  process,  <I>  is  a  normalized  correlation  matrix  with  the  Hermitian-Toeplitz  form 
where  each  diagonal  element  is  equal  to  one. 
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If  the  complex  amplitude,  a,  were  known  then  the  Neyman-Pearson  likelihood  ratio  (LR)  statistic 
is  given  by 


X  = 


(4) 


Define  the  function  FN  as 


(5) 


Then  it  is  seen  that  the  LR  detector  is  given  by 


X  =  F^(z  ~  aso)H  (z  ~  flSo»  (6) 

Fn( zh  4>  '  z) 

For  many  detection  problems,  the  complex  signal  amplitude  of  the  desired  signal  vector  is 
unknown  to  the  detector.  In  general  if  one  has  unknown  parameters  in  the  LR  detector,  a  scheme  known 
as  the  Generalized  Likelihood  Ratio  Test  (GLRT)  is  commonly  used  whereby  the  unknown  parameters 
are  found  using  their  maximum  likelihood  (ML)  estimates  under  each  hypothesis  [13].  These  estimates 
are  then  inserted  into  the  respective  pdfs  (under  H0  and  //,)  and  the  LR  is  formed.  We  now  derive  the 
GLRT  for  detection  of  a  signal  with  unknown  complex  amplitude  in  additive  Gaussian  mixture  noise. 

The  GLRT  for  our  problem  is  given  by 


Fn((z  -  <5s0)H  <f>-‘  (z  -  ds0» 
Fa(zh  <!>-'  z) 


(V) 


where 


a  =  arg  max  FN(( z  -  as0)H  4*'1  (z  -  as0)). 

a 


(8) 


The  maximum  of  FN  with  respect  to  a  can  be  found  by  setting  dFN/da  =  0  and  solving  for  a.  Now 


dFN  _  d(i  dFN  =  q  (9) 

da  da  dp 


where 


P  =  (z  -  as0)H  <I>-1  (z  -  as0). 


(10) 


Using  (5),  we  can  show 


3 


(11) 


dFN 

d($ 


\ 


1 


exp 


0 

T 


dP. 


From  (11)  it  is  seen  that  dFfJd(i  <  0.  Thus  all  solutions  for  a  in  Eq.  (9)  must  satisfy  dfi/da  =  0.  It 
is  straightforward  to  show  that  this  equation  has  only  one  solution  which  is  given  by 


s0H  <J>'z 
So” 


(12) 


Note  that  (12)  is  the  ML  estimate  of  the  complex  signal  amplitude  when  the  contaminating  noise  is  purely 
Gaussian.  Inserting  (12)  into  (7)  results  in  the  GLRT  for  Gaussian  mixtures: 


F^z"  z  (1  -  |p|2)) 
F^ z"  4»-'  z) 


where  |-|  denotes  magnitude, 


M2 


| So" 

(zH  <h-'  z)  (SoH  <h-'  s/ 


(13) 


(14) 


and  0  <  |p|  <1. 

B.  Problem  Statement 

Let  P  be  the  set  of  all  univariate  distribution  functions  that  have  support  fl  9[0,oo],  are  right- 

continuous  at  r  =  0,  and  equal  zero  at  r  =  0.  Let  G  and  G  be  non-negative,  right-continuous,  bounded 
monotonically  increasing  functions  (not  necessarily  distribution  functions)  on  [0,oo]  which  are  zero  at  r 
=  0.  Define  the  uncertainty  class  Pr  as 


yr  =  {Pr\G  <  Pt  <  G,  Pr  £  P}.  (I5) 

Often  one  implements  a  detection  scheme  based  on  a  representative  element  of  PT,  say  P,(l);  i.e. 

the  GLRT  given  in  the  previous  subsection  is  implemented  with  Pt  =  P,(l>.  Pr(l)  will  be  called  the 
detector  operating  point.  If  the  actual  multivariate  inputs  had  an  underlying  distribution  Pr<2)  £  PT,  we 
are  interested  in  formulating  the  detection  and  false  alarm  probability  performance  which  we  denote  by 
PD(Pr(1),  Prp),  7)  and  P^fP/0,  P,<2),  7),  respectively.  In  addition,  we  would  like  to  find  nontrivial  upper 
and  lower  bounds  on  these  performance  measures  such  that 


and 


Ld(G,G, 


T)  <  PD(Pr<u,P,<2>,7)  <  Ud(G,G,T) 


(16) 
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where  the  upper  and  lower  bounds  are  functions  of  just  G,  G,  and  the  detection  threshold,  T.  The 
detection  measures,  PD,  LD,  and  UD,  are  also  a  function  of  the  signal-to-noise  ratio.  In  the  following 
sections,  expressions  for  PD,  PF,  LD,  UD,  LF,  and  UF  are  derived. 

III.  PD  AND  PF  DERIVATION 

A.  Equivalent  GLRT 

In  this  subsection  we  present  a  simplified  equivalent  GLRT  which  will  allow  us  in  the  next  section 
to  find  explicit  results  for  the  GMD  detector’s  performance.  To  this  end,  it  is  known  [16]  that  if  4>  is 
non-singular,  then  we  can  matrix  transform  the  input  vector,  z,  by  a  non-singular  N  x  N matrix  A,  which 
has  the  properties  that  the  input  noise  vector  is  whitened,  each  input  element  has  noise  power  normalized 
to  one,  and 


As0  =  ((s0H  $->  s0)1/2,  0,0, . . .  ,0)T  *  s0.  <18) 

We  note  that  all  of  the  desired  signal  has  been  placed  into  the  first  element  by  the  matrix  transform. 
Thus  if  u  =  Az  where  u  =  («„  u2,  ,  %)T,  then  under  H0 


E  {uuH}  =  IN  (^9) 

where  IN  is  the  N  x  N  identity  matrix.  Transforming  the  input  by  a  non-singular  matrix  transform,  A, 
does  not  change  the  performance  measures  of  the  GLRT  that  were  discussed  in  the  preceding  section. 

If  we  use  the  vector,  u,  as  our  starting  point  in  determining  the  GLRT,  it  can  be  shown  that  the 
N-variate  pdf  of  u  under  the  two  hypothesis  is  given  by  (2)  and  (3)  by  setting  $  =  IN.  Furthermore,  the 
GLRT  is  given  by  (13)  and  (14)  with  4>  =  IN  as 


.  =  /y uH  u  (1  -  |p |2) 
Fn  (uh  u) 


where  FN  is  defined  by  (5)  and 


IpI2  = 


(uH  u)  (s0H  s0) 


N 


E  Kl2 


(20) 


(21) 


Substituting  (21)  into  (20)  results  in 


(22) 


E  kl2 

n -2 


E  k 


n  =  1 


Define  the  N-\ ariate  function,  A,  such  that  X  =  A(u).  A  functional  block  diagram  of  the 
equivalent  GLRT  is  shown  in  Fig.  1.  In  this  figure,  T is  the  detector  threshold. 

B.  Derivation 

In  this  section,  we  derive  closed-form  solutions  for  the  probabilities  of  detection,  PD  (Prw,  P}2), 
7),  and  false  alarm,  PF(Pfw,  P,a\  7),  associated  with  the  GLRT  discussed  in  the  previous  section.  After 
the  A  matrix  transformation,  the  desired  signal  is  contained  only  in  the  first  element  of  u  under  Hu  and 
is  equal  to  a(s0H  4>‘‘  s0)1/2.  We  set  a  =  a0  ^  where  </>  is  the  phase  and  a0  the  magnitude  of  a.  Note  that 
at  this  point  in  processing,  the  first  element  of  u  is  the  matched  filter  output  and  that  the  signal-to-noise 
(S/N)  associated  with  the  matched  filter  denoted  by  (S/N)^  can  be  expressed  as 


(S/NV 


«o2 


So  *- 


V 


(23) 


Set  a0  =  a0(s0H  #■'  s0)in.  We  now  analyze  the  statistically  equivalent  GLRT  presented  in  the  previous 
subsection.  Set 


(24) 


v  .  £  KP.  (25) 

2 

We  note  that  the  random  variables  (r.v.’s)  |«,  |2,  |«2|2, . . . ,  |ww|2  when  conditioned  on  t  are  statistically 
independent.  Thus  tv  and  v  when  conditioned  on  r  are  statistically  independent.  It  is  straightforward 
to  show  that  v  when  conditioned  on  r  under  H0  or  H]  has  a  2 (N  -  1)  Chi-square  pdf  [17]  with  a2  =  . 5 
and  is  given  by 


PWr(v|T)  =  — J -  vw-2  (26) 

V|T  7^-*  (7V-2)! 

where  p  |r  (-|x)  denotes  that  the  pdf  is  conditioned  on  t. 

Under  H0  and  conditioned  on  r,  u,  is  a  complex  circular  Gaussian  r.v.  with  variance  equal  to  t 
and  under  //,,  |m,|  has  a  non-central  Rayleigh  distribution.  In  fact,  it  is  straightforward  to  show 
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pw\Mt)  = 


1  e-w,T 
r 

_1  e  -(w  *  ti) 

T 


2  ajw 


,  w  >  0,  under  H0 


,  w  >  0,  under  H. 


(27) 


where  Ia  is  the  modified  O’th  order  Bessel  function  of  the  first  kind. 


The  decision  rule  in  terms  of  w  and  v  is 


F$\v) 


X  = 


H\ 

> 


F$\w  +  v) 


T. 


(28) 


where  the  superscript  on  the  FN  indicates  that  FN  is  evaluated  using  (5)  with  PT  =  PTW.  It  was  shown  in 
the  previous  section  that  FN  is  strictly  monotonically  decreasing.  As  a  result,  since  T  >  1,  and  because 
Fn  is  a  one-to-one  function,  an  equivalent  decision  rule  for  (28)  is 


w  <  F'n 
Hn 


<i)- 


v  =  T0(v,T) 


(29) 


where  FN'  is  the  inverse  function  of  FN.  Define  P(D  \  v,r)  to  be  the  probability  of  detection  when 
conditioned  on  v  and  r.  Thus 


P(D  |v,T) 


J. 

X 


e 


-(w 


dw 


( 


=  <? 


2 

2T0  (v,  T ) 

Nx’\ 

*  J 

(30) 


where  Q(-,  •)  is  Marcum’s  Q  function  [13].  Finally, 


PD  (P(T'\  Pr(2),  D  =  }n  }o"  P(D  |v,r)  p„|r(v|  r)  dvdP? . 


(31) 


The  false  alarm  probability  can  be  found  by  setting  a0  =  0  in  (30).  It  will  be  found  that 


PF  (Pr(1),  Pj2),T)  =  Jo  Jo"  pv|r(v|r)  exp 


F, S,"-1  (Fj}\v)/  T)  -  v 


dvdP; 


(2) 


(32) 
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IV.  PERFORMANCE  BOUNDS 


A.  Useful  Lemmas 

In  this  section,  we  prove  three  lemmas  that  will  be  used  to  derive  upper  and  lower  bounds  of 
performance  for  PD(P?\  P?\T)  and  P^P^,  Pf\T). 

Lemma  1 :  Let  fij)  be  a  differentiable  function  on  [0,oo]  with  J{0)  <  oo  and  A  oo )  =  o.  Define 


GJJ) 


G(t)  if  -  f(T)  >  0 
G(r)  otherwise 


(33) 


Gl(t)  = 


m 

G(t) 


if 


-f(T)  £  0 
otherwise . 


(34) 


where  '  denotes  differentiation  with  respect  to  t.  If  PT  E  !PT,  then 


fo“  At)  dGL  <  JQ  At)  dPr  <  j~  At)  dG, 


(35) 


Proof :  Integrating  by  parts,  we  can  show 


Jn  At)  dP,  =  jo°°  At)  dP,  =  -  /'(t)P,(t)  dr.  (36) 

Using  (33)  and  (34)  and  the  fact  that  G  <  Pt  <  G,  it  follows  that 

-  jo“  f(T)  Gl(t)  dT  <  -  Jo"  f(T)  Pr(T)  dT  <  -  |o”  f(T)  GJj)  dT  .  (37) 

Integrating  each  of  the  above  three  integrals  by  parts  results  in  (35).  □ 

Lemma  2:  Define 


m 


i 


or 


e-0,a  dG  + 


-L  e  -0,a  dG 

0/N 


(38) 


m 


r  fiw 

J  0 


1 


or 


^  ~P,Q 


dG 


+  f  "  _L  e  -*'«  dG 

J  0/n 


(39) 


Then 
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m  ^  w  ^  m- 


(40) 


Proof:  Integrating  by  parts,  we  can  show 


Set 


w 


f  ” 

d 

'  e -Vo' 

jo 

da 

aN 

PT(r)  dr  . 


(41) 


and 


-  (e  ~eia/aK)'  >  0  ('  denotes  d/da) 
otherwise 


(42) 


-  (e'^/ay  >  0 
otherwise. 


(43) 


It  is  straightforward  to  show  that  -{e  ~pialaN)'  >  0  iff  a  >  /3/A.  Also  from  (42)  and  (43),  it  follows 
that 


f  °° 

^3 

1 

1 _ 

'  e  -p>a' 

Jo 

i  da 

aN 

V(a)da  < 


d 

■  e  -m  • 

da 

aN  . 

Pr(T)dT  <  |o°° 


d 

'  e-**' 

da 

aN 

V(a)  da 


(44) 

Integrating  each  of  the  above  integrals  by  parts  and  using  the  definitions  given  by  (38)  and  (39),  Eq.  (40) 
results.  □ 

Using  the  definitions  of  FN,  F,  and  F  given  by  (5),  (38),  and  (39),  respectively,  it  can  be  shown 
that  Fn,  F,  and  F  are  strictly  monotonically  decreasing  functions  of /?  on  [0,oo].  Let  0^,  andfy. 
be  the  ranges  of  FN,  F,  and  F  respectively.  Each  of  these  ranges  includes  0  as  one  end  point  (/3  =  oo). 

The  other  end  point  of  the  ranges  may  or  may  not  be  infinite.  Now  FN\  F  \  and  F1  are  strictly 
monotonically  decreasing  functions.  We  define  on  [0,oo]  the  following  functions: 


<(X)  = 


F-'(X) 

0 


if 


X  e 

otherwise 


(45) 
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fo-'W  = 


F'(X) 

0 


if 


X  e  ftf 
otherwise 


(46) 


and 


^.o(X)  = 


FN\\)  if 
0 


X  6 

otherwise . 


Lemma  3:  For  T  >  1 


(F(v)/7)  <  F„'  (F„(v)/7)  <  F0'  (F(v)/T) . 


(47) 


(48) 


Proof:  For  F  >  1,  Fw  (v)/Fis  always  in  the  range  of  FN.  From  Lemma  2 ,  F  <  FN  <  F.  Using  this 
and  the  fact  that  F,  F,  and  FN  are  strictly  monotonically  decreasing,  it  can  be  shown  that 

=s  it 

Thus 


FjoiF^VD  ^  r'  (FA(v)/7)  >  F1  (F(v)/7)  (49) 


and 


F„,'o  (F,,(v)/7)  <  F0‘  (FAfv)/7)  <  F0‘  (F(v)/7).  (50) 

The  last  inequality  in  (49)  follows  because  F^'  is  monotonically  decreasing  and  FA(v)/F  <  F(v)/F.  A 

similar  argument  establishes  the  last  inequality  of  (50).  Since,  for  F  >  1,  FNl0  (FA(v)/F)  =  Fn'(FJ,v)/T) , 
and  the  lemma  follows.  □ 

B.  Derivation 

In  this  subsection,  we  derive  nontrivial  lower  and  upper  bounds  of  performance  for  P0(Pr(1),  Pr(J), 
7)  and  P^P/1’,  Pr(2),  7). 

Theorem  1: 


Ld(G ,  G,7)  <  PD{Pl'\P?\T)  <  U0(G,  G,7) 


(51) 


where 
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(52) 


Ld(G,  G,T)  =  Jo“  Hd(j, 7)  dGw 


GJj)  = 


Gud(t)  =  i 


,  G,T)  =  Jo”  Hd(t,T)  dGUD 

(53) 

JD  [Fi  (F(v)/T)  -  v,  T]  py|T(v|r)  dv 

(54) 

JD  ^  (W7)  -  v,  r]  py|T(v|r)  rfv 

(55) 

G(r)  if  -  dHJdr  >  0 

(56) 

G(r)  otherwise 

G{t)  if  -  dHJdr  >  0 

(57) 

GJr)  otherwise 

/ 

_ 

\ 

2 

2 1 

t’  \ 

)(2) 


and 

JD(t,  t)  =  /*“  -  e^w  *“o)h  I0(2a0\fw  /t)dw  =  <? 

Prao/:  From  the  results  of  Section  III,  it  can  be  demonstrated  that 

PD{P?\  P?\T)  =  Jfi  Hd\t,T)  dP, 

where 

Hd\t,T)  =  [  ”  f  “  -  e'<H'*a”)/T  I  (2a0\fw It)  pv[t(v\t )  dwdv, 

J  0  J  T0(v,T)  T 

pvjT  is  given  by  (26),  and  T0(v,T)  is  given  by  (29).  Furthermore, 

Hd\t,T)  =  J"  JD  [(FNl(FJy)/T)  -  v,  r]  p,|T(v|r)  dv. 


(58) 


(59) 


(60) 


(61) 
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Using  Lemma  3,  since  JD  is  monotonically  decreasing  with  respect  to  its  first  argument,  we  can 


show 


-  V,  t)  <  JD  [Fn'(F^v)/D  -  v,  r] 
ss  (F(v)/T)  -  v,  7], 


(62) 


Since  pv|t(v|t)  0,  it  follows  that 

Hd(t,T)  <  Hq\t,T)  <  Hd(t,D  (63) 

and 

|o“  Hd(t, 7)  dPf]  <  |o  hS\t,T)  dP r<2)  <  |o“  W0(r,r)  rf/>,(2).  (64) 

It  is  straightforward  to  show  that  H_d  and  HD<  oo  for  7G[0,  o°],^D(°° ,7)  =  7/D(°°,7)  =  0,  and  // 
and  //D  are  differentiable  on  [0,  oo].  Thus  invoking  the  appropriate  inequalities  of  Lemma  1: 


J0‘  UJr.T)  dGw  <  jj  H„<r,7)  dP ®  (65) 

and 

|o”  Hd(t>T)  dPl2)  <  |o”  Hd(t,T)  dGVD  (66) 

where  G^,  and  G^  are  defined  by  (56)  and  (57),  respectively.  The  results  of  (64)-(66)  prove  the 
theorem.  □ 

In  similar  fashion  we  can  prove 
Theorem  2: 


L^G,  G,T)  <  P^\  P?\T)  <  UfiG,  G,T) 

where 


Lf(G,  G,T)  = 

J  ’  « ('.n 

G/G,  G,7)  = 

f0‘  "An 

(67) 


(68) 


(69) 
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KF(r,T) 


f  00 

F0'(F(v)/T)-v 

Jo  eXP  - 

r 

Pv|r(VlT)  dV 


(70) 


and 


H^r,T) 


r  00 

F^KF(v)/T)-v 

Jo  eXP  ■ 

T 

pM7)  dv 


(71) 


GJ7) 


G(r) 

_G{t) 


if 


-dH  /dr  >  0 

— f 

otherwise 


(72) 


Gurkr) 


Gir) 

G(r) 


if 


-dHJdi  >  0 
otherwise . 


(73) 


C.  Performance  Spread 

For  a  given  detector  operating  point  and  threshold,  it  is  desirable  to  quantify  the  performance 
measure  spread  over  the  input  uncertainty  class.  Thus  we  define 


Sd(P(o\T)  = 


max  |  PD(PT 

\  P?}  €  7 


(1) 


P?\T)  -  PD(Pl'\  P?\T)\ 


>(3) 


(74) 


and 


S^P(f'\T) 


max 

pQ)  pO)  £  ep 


\P/Pr'\  P?\T)  ~  PF(P?\  P?\T)\ 


>0» 


(75) 


We  desire  to  upper  bound  the  maximum  performance  measure  spread.  To  this  end,  define 


Sd(T)  =  max  •Sd(/\<1),7) 

/>f<0  G  y 


(76) 


S^T)  =  max  S^t'\T). 

Ft(,)  €  7 

Using  Theorems  1  and  2  it  is  trivial  to  show  that 


Sd(T)  <  Ud(G,  G,T)  -  Ld(G,  g,d 


(77) 


(78) 
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(79) 


SJJ)  ^  UJG,  G  T)  -  LfG,  G,T). 


The  following  theorem  provides  a  looser  set  of  upper  bounds  on  SD  and  Sr 

Theorem  3:  Let  G  and  G  admit  density  functions  g  and  g,  almost  everywhere  with  respect  to  measure, 
r  (a.e.  [t])  respectively  where  g  =  dG/dr,  g  =  dG/dr.  Then 

Sd(T)  <  l~  Hd(t,T)  \g  -  g\dT  (80) 


SJT)  ^  }o“  \g  -  g\dr. 


(81) 


Proof :  The  proofs  of  (80)  and  (81)  are  similar  so  that  we  only  present  the  proof  of  (80).  From  (51)-(53), 
it  follows  that 


Sd(P?\T)  £  jo“  h2\t, T)  d{Guu  -  Gu). 


(82) 


If  G  and  G  admit  densities,  then  Gw  and  GUD  defined  by  (56)  and  (57)  admit  densities  a.e.  [r].  Set 
8ld  =  dGJdr,  gUD  =  dGJdr .  In  addition,  \gUD  -  gw \  =  |g  -  g\.  Using  (82) 

W". T)  <  |  *  (gm  -  w  d7 

£  J  *  Hd(t,T)  \g„  -  gu>\  dr  -  Jo"  Hd(t,7)  |I  -  s|  dr.  (83) 

Since  the  expression  on  the  right  above  is  independent  of  Pt°\  (80)  follows.  □ 

The  bounds  given  above  might  be  utilized  as  follows.  First,  find  the  threshold  T,  such  that  a  given 

desired  probability  of  false  alarm  is  equal  to  Gf(G,  G,  T).  Hence  we  can  guarantee  the  false  alarm 
performance  over  the  uncertainty  class  for  any  detector  operating  point  that  is  chosen  from  the  uncertainty 

class.  Next,  find  UD(G,  G,  T)  and  LD(G,  G,  T)  and  ascertain  1)  is  LD  sufficient  to  guarantee  a  given 
level  of  detection  performance' and  2)  what  is  the  potential  spread  (Use  (78)  or  (80))  of  detection 
performance? 

V.  SUMMARY 

Detection  performance  in  terms  of  probabilities  of  detection  and  false  alarm  (Pn  and  PF)  was 
derived  for  a  desired  signal  embedded  in  correlated  multivariate  noise  from  a  Gaussian  mixture 
uncertainty  class.  This  uncertainty  class  is  defined  using  upper  and  lower  bounding  functions  on  the 
univariate  Gaussian  mixing  distribution  function.  For  a  given  detector  operating  point,  expressions  for 
PD  and  PF  were  derived.  In  addition,  nontrivial  lower  and  upper  bounds  of  performance  were  derived 
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which  are  functions  of  the  bounding  functions  of  the  uncertainty  class,  the  detection  threshold,  and  the 

matched  filter  signal-to-noise  power  ratio,  (S/N)^,  (only  for  PD). 
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Fig.  1  —  Equivalent  detector 
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